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Ax  —  1,  «  =  1.  D  =  0.20, 
k  infinite 


ABSTRACT 

Some  of  the  observed  complexities  of  drying  rate  curves  have  been 
ascribed  to  a  variability  of  the  diffusivity  of  the  substance  with  changing 
moisture  content.  Quantitative  evaluation  of  the  suggestion  has  been 
difficult  because  the  generalized  diffusion  equation  is  not  integrable. 
Approximate  numerical  procedures  are  applied  in  the  present  paper  to 
explore  the  effect  of  a  variable  diffusivity  in  the  low-moisture  range  upon 
the  shape  of  the  rate  curve  for  the  infinite  slab,  the  long  cylinder,  and 
the  sphere.  Either  vapor  pressure  or  concentration  may  be  used  as  the 
driving  potential.  The  effects  of  changing  the  thickness  of  the  body  and 
changing  the  relative  importance  of  the  surface  mass  transfer  coefficient 
are  investigated.  Typical  experimental  rate  curves,  determined  during 
the  dehydration  of  cut  vegetables,  exhibit  curvature  of  the  type  predicted 
by  such  calculations,  and  to  a  degree  that  might  well  be  expected  on 
the  basis  of  independent  measurements  of  diffusivity  in  hydrophilic  solids. 

THE  phenomena  occurring  during 
the  drying  of  a  moist  solid,  ever 
since  the  classical  analyses  of  Lewis, 
Sherwood,  and  others,  (21)  (22)  (26) 
(27)  (31)  have  been  divided  into  two 
or  more  distinct  phases.  The  earlier 
of  these,  characterized  by  phenomena 
similar  to  those  of  an  evaporating 
free  liquid  surface,  and  known  as  the 
constant-rate  phase,  has  presented  few 
difficulties.  The  one  or  more  later 
phases  of  the  process,  characterized 
by  a  steadily  decreasing  drying  rate, 
exhibits  complexities  which  have  not 
yet  been  fully  unraveled.  Various 
physical  mechanisms  have  been  called 
upon  to  suggest  explanations  for  the 
phenomena.  Sherwood  (31)  and 
Newman  (27)  employed  the  diffusion 
theory  first  worked  out  for  the  un- 
steady-state conduction  of  heat,  con- 
sidering moisture  concentration  gra- 
dient as  being  analogous  to  tempera- 
ture gradient.  Ceaglske  and  Hougen 
(P)  and  Fisher  (13)  emphasized  the 
mechanics  of  flow  through  capillary 
pore  spaces.  Macey  (23a)  deter- 
mined the  pressure  necessary  to  re- 
move water  from  moist  clays,  and 
identified  the  diffusional  driving 
force  with  this  pressure.  McCready 
and  McCabe  (26),  and  Bateman, 
Hohf,  and  Stamm  (5)  postulated  that 
vapor  diffusion  occurs  within  voids 
and  considered  that  two  or  more  of 
these  mechanisms  were  effective  si- 
multaneously. 

During  the  past  five  years  workers 
in  this  laboratory  have  made  a  com- 
prehensive investigation  of  the  drying 
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conditions,  such  as  the  use  of  com- 
mercial types  of  drying  trays  and 
heavy  tray  loadings,  which  make 
theoretical  analysis  of  the  results  dif- 
ficult. The  methods  used  and  the  re- 
sults of  the  work  have  been  reported 
by  Brown  and  Kilpatrick  (6)  and  in 
a  series  of  mimeographed  circulars. 
(36)  .  . 

Government  purchase  specifications 
for  dehydrated  vegetables  required 
that  the  moisture  content  be  very  low 
(4  to  6  per  cent,  depending  upon  the 
product).  The  drying  rates  in  the 
lower  part  of  the  falling-rate  phase 


rates  of  the  five  or  six  vegetables 
which  were  most  important  in  the 
wartime  dehydration  program.  The 
primary  purpose  of  this  work  was  to 
provide  the  data  needed  for  rational 
design  of  vegetable  dehydrators.  Con- 
sequently most  of  it  was  done  under 


therefore  took  on  great  practical  sig- 
nificance. Experimental  results  show, 
as  pointed  out  by  Van  Arsdel,  (37) 
that  as  much  time  may  be  required  to 
reduce  the  moisture  content  of  a  vege- 
table from  7  per  cent  to  4  per 
cent  as  is  required  to  remove  the 


FIG.  2.  MOISTURE 
GRADIENTS,  CONSTANT 
PERMEABILITY 

Slab,  half -thickness  4Vj, 
Az=l,  t=l,  P  —  4.40 
X  10-«,  k  infinite.  Vapor 
pressure  isotherm  as  in 
Fig.  3 
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PIG.   3.     TYPICAL  VAPOR-PRESSURE 
ISOTHERM 

Dehydrated  carrots,  140  Degrees  P. 

hundred-fold  larger  amount  of  mois- 
ture from  the  starting  point  down  to 
7  per  cent.  Since  the  size  of  the  drier 
required  is  directly  proportional  to  the 
drying  time  for  a  specified  rate  of 
throughput,  the  design  significance  of 
the  final  phase  is  obvious. 

During  the  empirical  analyses  of 
experimental  drying  data  in  this  lab- 
oratory it  was  noted  that  drying 
times  of  vegetable  pieces  could  be 
correlated  satisfactorily  by  means  of 
nomographs  provided  the  data  were 
regarded  as  falling  into  two  distinct 
phases,  the  dividing  line  being  in  the 
range  between  10  and  20  per  cent 
moisture,  (d)  In  the  low-moisture 
phase  the  characteristics  of  the  data 
corresponded  with  those  of  a  system 
in  which  the  rate  of  internal  diffusion 
is  controlling;  that  is,  neither  the  air 
velocity,  the  type  of  tray,  nor  the  orig- 
inal depth  of  tray  loading  had  any 
further  appreciable  effect.  The  data 
should  be  amenable  therefore  to  theo- 
retical analysis  if  the  mechanism  of 
moisture  transfer  in  nearly  dry  mate- 
rials of  this  character  is  similar  to  any 
of  those  postulated  by  previous  inves- 
tigators. The  present  study  is  directed 
toward  the  moisture  diffusion  hypoth- 
esis of  Sherwood  and  Newman,  and 
certain  modifications  of  it,  rather  than 
the  pore  flow  theory  of  Ceaglske  and 
Hougen,  because  the  materials  in 
question  resemble  solid  gels  in  struc- 
ture. 

In  this  classical  diffusion  theory  the 
potential  producing  movement  of 
moisture  within  the  solid  piece  is  iden- 
tified with  a  concentration  difference. 
The  differential  equation  expressing 
the  change  in  concentration  at  any 
point  with  respect  to  time  is  as  fol- 
lows for  one-dimensional  diffusion 
(the  infinite  slab)  : 


where  D  is  the  diffusivity  of  moisture 
in  the  material.  Numerous  writers 
have  examined  the  integrals  of  this 


equation  for  various  important  boun- 
dary conditions,  (3)  (7)  (15)  (16) 
(21)  (27)  (31)  (32)  and  the  pre- 
dicted drying  rates  have  been  shown 
to  agree  reasonably  well  with  experi- 
mental data  on  the  drying  of  some  in- 
dustrial solids.  However,  when  the 
experimental  rate  curves  are  ex- 
amined in  detail,  divergence  from  the 
theoretical  shape  is  apparent  in  nearly 
every  case.  The  divergence  is  most 
apparent  in  the  range  of  very  low- 
moisture  contents,  where  the  solu- 
tions of  Equation  ( 1 )  predict  a  direct 
proportionality  between  drying  rate 
and  free-moisture  content.  Sher- 
wood, (31)  Bateman,  Hohf,  and 
Stamm,  (5)  Hougen,  McCauley,  and 
Marshall,  (16)  and  King  (20)  con- 
clude that  diffusivity  does  not  remain 
constant,  but  is  a  function  of  moisture 
content.  The  mathematical  complex- 
ity of  the  resulting  case  is  such  that 
no  very  useful  analytical  solutions  for 
the  diffusion  equation  have  been 
found,  and  the  observation  has  re- 
mained a  purely  qualitative  one. 

There  is  now  ample  evidence  that 
the  diffusivity  of  moisture  in  hydro- 
philic  substances  indeed  is  highly  de- 
pendent upon  the  moisture  content. 
The  results  of  Doty,  Aiken,  and  Mark 
(10)  indicate  that  the  diffusivity  of 
moisture  through  regenerated  cellu- 
lose changes  by  at  least  forty-fold 
within  a  comparatively  narrow  range 
of  moisture  contents.  King,  (20) 
studying  the  moisture  relations  of 
keratin  proteins,  such  as  horn,  founc 
more  than  a  hundred-fold  change. 
The  failure  of  experimental  drying 
curves  to  follow  the  course  predicted 
by  the  simple  diffusion  theory  is 
scarcely  surprising. 

Approximate  numerical  solution  of 
the  diffusion  equation  offers  no  great 
difficulty,  and  it  possesses  the  advan- 
tage that  some  of  the  simplifying 
assumptions  that  are  necessary  to  the 
analytical  integration  of  Equation 
(1)  may  be  dropped  if  desired.  The 
two  most  important  of  these  assump- 
tions, from  our  viewpoint,  are  the  fol- 
lowing : 

1.  The  potential  producing  diffu- 
sion is  a  concentration  difference. 

2.  The  diffusivity  is  independent  of 
the  moisture  content. 

Some  other  assumptions  that  must 
be  introduced  into  the  analytical  treat- 
ment become  unnecessary  in  the 
numerical  method.  Such  is  the  case, 
for  example,  with  respect  to  the 
shrinkage  in  dimensions  as  moisture 
is  lost,  to  the  existence  of  tempera- 
ture gradients  within  the  piece,  and 
to  nonhomogeneity  of  the  material 
and  variation  of  diffusivity  with  di- 
rection in  the  piece.  The  limiting  fac- 
tor, in  most  practical  cases,  is  likely 
to  be  the  difficulty  of  experimental 


determination  of  the  necessary  data. 
These  further  complexities  are  not 
considered  in  the  present  study. 
Shrinkage  and  the  existence  of  tem- 
perature gradients  may  be  neglected 
in  fact  without  appreciable  error  if 
attention  is  directed  solely  to  the  re- 
gion of  low-moisture  content. 

The  Driving  Force  for  Diffusion 

Equation  (1)  is  derived  from  an 
expression  analogous  to  Newton's  law 
of  cooling  and  to  Ohm's  law.  The 
common  characteristic  of  these  two 
laws  is  that  they  equate  a  rate  of  flow 
to  a  certain  potential  difference  di- 
vided by  a  resistance  (or  multiplied 
by  a  conductance).  The  potential  dif- 
ference is  pictured  as  an  unbalanced 
force,  producing  motion  against  some- 
thing comparable  to  friction. 

Several  characteristics  of  moist 
bodies,  besides  concentration  differ- 
ences, might  be  pictured  in  the  role  of 
diffusion-producing  potentials.  An 
obvious  one  is  the  activity  of  the 
moisture  at  any  given  point,  as  meas- 
ured by  its  equilibrium  vapor  pressure 
there.  The  Ohm's  law  type  of  expres- 
sion still  can  be  used,  provided  the 
conductance  term  is  written  also  in 
appropriate  new  units.  It  is  now  cus- 
tomary in  diffusion  studies  to  name 
this  conductance  variable  the  per- 
meability when  the  potential  is  vapor 
pressure,  while  diffusivity  is  the  cor- 
responding term  when  the  potential  is 
concentration.  Then  (see  McCready 
and  McCabe  (26) )  the  general  equa- 
tion (written  to  allow  for  a  variable 
permeability)  will  become : 

A£  =  A.(pAPl)  (2) 

Or,  in  the  notation  more  commonly 
used  in  drying  studies, 

iZl  =  A    U  p  dL\ 

30      s      dx\  dx) 

(3) 

It  must  be  emphasized  that  the 
choice  of  potentials  and  conductances 
is,  in  the  absence  of  a  temperature 
gradient,  only  a  matter  of  conven- 
ience; if  the  proper  units  are  con- 
sistently used  the  solutions  of  the 
basic  equation  are  identical.  This  is 
not  to  say,  however,  that  one  choice 
may  not  be  better  than  another  theo- 
retically. It  is  to  be  expected  that 
some  one  of  the  hypothetical  mechan- 
isms will  be  the  most  nearly  in  accord 
with  other  physical  knowledge  about 
the  substance.  Three  different  lines 
of  investigation  rray  be  pertinent  to 
this  theoretical  question. 

Consider  first  the  temperature  co- 
efficient of  diffusion  rate.  If  concen- 
tration is  considered  as  the  potential 
then  the  entire  effect  of  a  change  in 
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temperature  of  the  system  must  be 
ascribed  to  a  change  in  the.diffusivity 
of  the  solid.  If  vapor  pressure  is  the 
potential,  on  the  other  hand,  some  of 
the  change  in  rate  of  transfer  must 
be  ascribed  to  a  change  in  the  poten- 
tial difference  and  much  less  to  a 
change  in  the  properties  of  the  solid. 
Measurements  of  the  effect  of  tem- 
perature on  swelling,  molecular  spac- 
ing, or  other  properties  may  support 
one  of  these  conclusions  against  the 
other. 

Or  again,  consider  the  variability 
of  diffusivity  (and  permeability)  with 
moisture  content.  In  general,  the 
water  vapor  pressure  of  hydrophilic 
substances  does  not  obey  Henry's  law. 
The  two  constants,  D  and  P,  are  re- 
lated to  each  other  through  the  actual 
vapor-pressure  isotherm  of  the  sub- 
stance,* and  in  all  the  common  hydro- 
philic substances  that  isotherm  is 
strongly  curved,  generally  sigmoid  in 
shape.  The  consequence  of  this  cur- 
vature is  that  if  either  D  or  P  is  inde- 
pendent of  moisture  content,  the  other 
must  vary  with  moisture  content ;  or, 
more  generally,  the  variations  of  D 
and  P  with  moisture  content  will  fol- 
low quite  different  courses.  Once 
again,  measurements  of  other  physi- 
cal properties  of  the  substance  may 
accord  better  with  one  of  these  pic- 
tures than  with  the  other. 

Diffusion  in  nonisothermal  solids 
has  received  only  the  most  superficial 
study.  The  few  observations  (1) 
(3)  (23)  are  only  semiquantitative, 
and  the  conclusions  drawn  are,  in 
part,  contradictory.  It  does  appear  to 
be  certain,  however,  that  moisture 
transfer  occurs  within  a  uniformly 
moist  body  under  the  action  of  an  im- 
posed temperature  gradient.  In  this 
case,  surely,  concentration  difference 
may  not  be  regarded  as  the  potential 
term,  for  there  is  no  concentration 
difference.  Apparently,  the  converse 
case  never  has  been  investigated ;  con- 
sider, that  is,  whether  steady  flow  of 
moisture  through  a  diffusion  barrier 
would  occur  if  a  temperature  gradient 
were  maintained  in  the  barrier  but  the 
vapor  pressure  of  moisture  were 
maintained  the  same  on  both  faces  of 
the  sheet.  In  the  latter  case  there 
would  be  no  vapor-pressure  gradient, 

*  If  the  units  used  are  consistent,  the  con- 
stants are  related  as  follows : 


D  = 


A  consistent  set  of  units  would  be,  for  ex- 
ample, for  D :  sq.cm.  per  sec. ;  for  5 :  g.  per 
cu.cm. ;  for  p :  mm.  Hg. ;  for  T:  dimension- 
less  ;  for  P :  g.  per  sec.  through  1  sq.cm.  un- 
der a  gradient  of  1  mm.  Hg.  per  cm.  thick- 
ness (or  (g.)/(sec.)  (sq.cm.)  (mm.Hg./ 
cm.)).  Since  the  slope,  dp/dT,  of  the  va- 
por-pressure isotherm  changes  with  mois- 
ture content,  the  relation  between  D  and  P 
is  likewise  variable. 


FIG.    4.      DRYING    RATES,  CON- 
STANT   DIFFUSIVITY    AND  CON- 
STANT PERMEABILITY 

Slab,  half  thickness  4V4,  Az  =:  1,  «  =  1. 
Curve  A,  D  =  0.20,  k  infinite.  Curve 
B,  P  =  4.40  X  10-«,  k  infinite,  vapor 
pressure  as  in  Fig.  3 


but  there  would  be  in  general,  a  con- 
centration gradient  across  the  sheet. 
Certain  of  the  lines  of  investigation 
referred  to  in  this  and  the  preceding 
paragraphs  are  being  pursued  by 
other  workers  in  this  laboratory. 

If  the  mechanism  of  diffusion  is 
pictured  in  the  manner  suggested  by 
Barrer   (3)    (4)   and   Stearn  and 
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FIG.    5.     MOISTURE    CONTENT  AT 
CENTER-LINE 
Conditions  as  in  Fig.  4.   Curve  A,  dif- 
fusivity constant.   Curve  B,  permeabil- 
ity constant 


Eyring  (33)  and  further  developed 
by  Cassie  (8)  and  King  (20),  the 
form  of  the  vapor-pressure  isotherm 
should  be  correlated  with  that  of  the 
permeability  vapor-pressure  relation. 
This  approach,  which  is  that  of  sta- 
tistical mechanics,  postulates  that 
water  molecules  which  have  absorbed 
more  than  a  certain  minimum  activa- 
tion energy  are  able  to  move  across 
energy  barriers  characteristic  of  the 
molecular  structure  of  the  solid.  Ap- 
plication of  the  Clausius-Clapeyron 
relation  to  the  vapor-pressure  iso- 
therms of  various  organic  materials 
shows,  however,  that  the  binding  en- 
ergy for  water  increases  substantially 
and  rapidly  at  low-moisture  contents  ; 
the  calculated  ratio  of  the  latent  heat 
of  vaporization  of  water  from  dehy- 
drated carrots  to  the  latent  heat  of 
vaporization  of  free  water  at  the  same 
temperature  rises  from  very  nearly 
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unity  at  8  per  cent  moisture  to  about 
1.60  at  2  per  cent  moisture.  The 
height  of  the  postulated  energy  bar- 
riers evidently  should  be  substantially 
greater  at  the  lower-moisture  level, 
and  the  permeability  therefore  would 
be  lower. 

Materials,  such  as  dehydrated  vege- 
tables, are  highly  complex,  multicom- 
ponent  systems.  Their  diffusion  char- 
acteristics must  be  expected  to  be  a 
composite  of  the  characteristics  of 
crystalline  or  glassy  sugars,  gelled 
starch,  coagulated  proteins,  cellulosic 
cell  walls,  pectins,  lipids,  and  other 
constituents.  Even  if  each  of  these 
constituents  should  exhibit  a  rela- 
tively sharp  break  in  properties  at 
some  level  of  moisture  content,  the 
combined  effect  of  all  of  them  would 
be  to  produce  a  more  gradual  change. 
To  the  extent  that  the  individual  con- 
stituents are  segregated,  rather  than 
in  solid  solution  in  one  another,  va- 
rious levels  of  moisture  content  may 
exist  side  by  side  in  equilibrium  be- 
cause of  the  different  activity  of  mois- 
ture sorbed  in  the  different  constit- 
uents. 

Approximate  Methods  of  Integration 

Approximate  methods  for  the  solu- 
tion of  the  diffusion  equation  have  re- 
cently been  the  subject  of  publications 
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FIG.   6.     WATER  VAPOR  PERMEABILITY 
OF   CELLOPHANE,   25   DEGREES  C. 

Deduced  from  data  of  Doty,  Aiken,  and 
Mark,(/0)  Permeability  in  cc.  water 
vapor  N.T.P./sec.  through  1  cm.2  for 
gradient    1    cm.    Hg/mm.  thickness 
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FIG.  7.    MOISTURE  DIFFUSIVITY  IN 
CELLOPHANE,  25  DEGREES  0. 

Deduced  from.  (10)  and  (17)  Diffusiv- 
,      ity  In  cm.Vsec. 

by  Dusinberre,  (11)  Paschkis  and 
Baker,  (2P)'  and  others.  (12)  (32) 
The  graphical  method  described  by 
Schmidt  (30)  was  shown  by  Dusin- 
berre to  constitute  a  special  case,  and 
to  possess  decided  limitations.  The 
modification  of  Schmidt's  method  to 
take  account  of  a  variable  diffusivity, 
referred  to  briefly  by  Sherwood  and 
Reed,  (32)  appears  from  unpublished 
studies  made  by  A.  H.  Brown,  of  this 
laboratory,  to  be  of  doubtful  general 
validity,  and  at  any  rate  to  be  inap- 
plicable if  the  diffusivity  varies  sub- 
stantially. A  semigraphical  method 
proposed  by  Dusinberre  to  handle  the 
variable  diffusivity  case  likewise  has 
been  found  unsatisfactory  if  the  dif- 
fusivity varies  widely  and  abruptly. 

Numerical  methods  of  solution  are 
based  upon  mass  balances  applied  to 
small,  but  finite,  elements  of  the  solid 
for  small,  but  finite,  increments  of 
time.  As  in  Dusinberre's  analysis,  we 
make  the  following  assumptions  : 

1.  A  time  interval  can  be  chosen 
sufficiently  small  that  there  is  a  negli- 
gible error  in  using  the  initial  poten- 
tial gradient  to  compute  the  flow  of 
moisture  during  this  interval.  (Note 
that  even  though  the  potentials  may 
be  changing  rapidly  in  two  adjoining 
elements,  the  difference  between  the 


two  potentials,  which  determines  the 
gradient,  changes  very  much  more 
slowly.) 

2.  The  time  interval  is  sufficiently 
small  that  there  is  negligible  error  in 
neglecting  the  effect,  during  that  in- 
terval, of  any  elements  beyond  the 
immediately  adjacent  ones. 

3.  The  dimensions  of  the  volume 
element  are  sufficiently  small  that 
there  is  negligible  error  in  using  the 
potential  at  its  center  point  to  com- 
pute the  moisture  content  of  the  ele- 
ment, or  vice  versa. 

Following  Dusinberre's  suggestion, 
we  find  it  convenient  to  make  the 
outermost  volume  elements  only  half 
as  thick  as  the  interior  ones,  so  that 
the  "center  points"  of  these  outermost 
elements  will  fall  on  the  surface  of 
the  body.  Film  coefficients  of  mass 
transfer  to  the  air  phase  are  then 
readily  introduced  into  the  calcula- 
tion. 

The  Working  Equations 

Consider  one-dimensional  diffusion 
in  an  infinite  slab.  Then,  using  vapor 
pressure  as  a  typical  potential  term, 
the  change  in  moisture  content,  AT,  in 
the  nth  element  from  the  surface  dur- 
ing the  time  interval  A0  will  be: 

AT„  =  [Pn+i,n(pn+i-  P*) 
AQ 


+  Pn-l,n(Pn-l  5(  A*)  2  ' 

(4) 

Since  the  permeability  may  undergo  a 
radical  change  within  the  range  of 
/>„_!  to. pn+i,  the  proper  evaluation  of 
the  terms  within  the  bracket  is  cru- 
cial. In  general, 

JfPm 
Pdp. 
Pn 

(5) 

The  permeability  will  be  determined 
usually  as  an  empirical  function  of 
moisture  content,  but  equally  well  may 


be  plotted  or  tabulated  as  a  function 
of  vapor  pressure  by  use  of  the  known 
moisture  content-vapor  pressure  rela- 
tionship. If,  Equation  (5)  is  put  into 
the  form, 


Table  1. 

Mean  moisture  content  of  infinite  slab.  Initial  moisture  content  uniform,  T  =  0.2000. 
Equilibrium  moisture  content,  T.  =  0.0200.  Flow  proportional  to  concentration  gradient. 
Diffusivity  constant,  D  =  0.20.  Surface  film  coefficient,  k  =  oo.  s  =  1,  A*  =  1.  Half- 
thickness  of  slab  divided  for  approximate  calculation  into  five  elements. 

Time 

Approximate 
Exact 
Error 


IS  30 
0.1221  0.0904 
0.1216  0.0903 

+.0005  +.0001 


50 
0.0632 
0.0632 
0 


60 
0.0538 
0.0538 
0 


70  80  85 

0.0464  0.0406  0.0382 

0.0465  0.0407  0.0.184 

-.0001  -.0001  -.0002 


Table  2. 


Mean  moisture  content  of  infinite  slab.  Initial  moisture  content  uniform,  T  =  0.2000. 
Equilibrium  moisture  content,  T.  =  0.0200.  Flow  proportional  to  concentration  gradient. 
Diffusivity  constant,  D  =  0.20.  Surface  film  coefficient,  k  —  0.10.  s  —  1,  Ax  =  1.  Half- 
thickness  of  slab  divided  for  approximate  calculation  into  5  elements. 


Time  30 
App.  0.1387 

Exact  0.1402 

Error  —.0015 


50 
0.1127 
0.1142 
—.0015 


75 
0.0681 
0.0891 
—.0010 


100 
0.0700 
0.0710 
-.0010 


120  140 
0.0590  0.0504 
0.0597  0.0510 

-.0007  —0.0006 


160 
0.0438 
0.0440 
—.0002 


180 
0.0385 
0.0388 
-.0003 


where  p0  is  any  convenient  base  value 
of  the  vapor  pressure,  the  value  of 
the  graphically  or  numerically  deter- 
mined integral  may  be  plotted  against 
vapor  pressure  and  the  terms  in  Equa- 
tion (4)  then  may  be  computed  read- 

The  change  in  moisture  content  in 
the  surface  element  may  be  similarly 

shown  to  be : 


ATi  =  2[P21(p2-pt) 

+  kAx(P*-p1)Y  ' 


s(Ax)1 


(7) 


Thus,  if  the  vapor  pressures  at  the 
reference  points  throughout  the  slab 
at  a  given  instant  of  time  are  known, 
the  change  of  moisture  content  at 
each  of  these  points  in  a  time  interval 
A6  can  be  calculated.  Two  different 
procedures  may  be  used  for  comput- 
ing successive  states  of  the  system: 
(a)  The  changes  in  T  may  be  fol- 
lowed from  step  to  step,  and  when 
the  new  set  of  values  of  T  has  been 
established,  corresponding  values  of 
p  may  be  read  from  the  vapor-pres- 
sure isotherm;  or  (b),  only  the 
changes  in  p  may  be  followed  from 
step  to  step;  a  curve  giving  dp/dT 
may  be  derived  from  the  isotherm; 
then : 


Ap  =  ^rAT. 


(8) 


The  second  procedure  is  generally 
preferable.  Values  of  T  throughout 
the  slab  do  not  need  to  be  evaluated 
at  every  step,  but  can  be  read  from 
the  isotherm  at  the  time  intervals  for 
which  they  are  wanted. 

Equation  (8)  is  relatively  precise 
for  small  values  of  AT,  but  a  second- 
order  correction  may  be  applied  with- 
out difficulty  if  desired;  instead  of 
reading  dp/dT  for  the  value  of  p  at 
the  beginning  of  the  time  interval,  it 
can  be  read  at  an  estimated  mean 
value  of  p  for  the  interval,  calculated 
mentally  by  subtracting  one-half  of 
the  value  of  Ap  found  for  the  pre- 
ceding time  interval. 

If  an  unduly  large  value  of  Ad  is 
used  in  the  stepwise  calculation,  or  an 
unduly  small  value  of  Ax,  the  calcu- 
lations will  exhibit  an  instability 
which  eventually  leads  to  absurd  r 
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suits.  It  is  as  though  a  resonant  oscil- 
lation, with  rapidly  increasing  ampli- 
tude, were  set  up  in  the  system. 
Dusinberre  (11)  proposes  the  use 
of  a  modulus  which  will  assure  a 
proper  choice  of  time  and  space  in- 
crements. In  general,  the  time  incre- 
ment can  be  increased  gradually  as 
the  system  approaches  equilibrium 
conditions.  The  resulting  saving  in 
labor  makes  it  worth  while  to  try 
this  expedient  at  intervals.  If  the 
oscillation  effect  then  begins  to  ap- 
pear, the  computer  can  return  to  the 
shorter  time  interval  and  the  diverg- 
ence will  be  damped  out  quickly. 

Lederer  (21)  and  Newman  (27) 
have  deduced  a  theorem  which,  at 
least  in  the  classical  diffusion  theory, 
makes  it  possible  to  pass  readily  from 
a  solution  for  a  one-dimensional  sys- 
tem (infinite  slab)  to  solutions  for 
corresponding  two-dimensional  (in- 
finitely long  rectangular  strip)  and 
three-dimensional  (rectangular  block) 
systems.  Unfortunately,  the  demon- 
stration appears  to  be  valid  only  for 
a  constant  diffusivity.  The  approxi- 
mate numerical  method  becomes  ex- 
cessively laborious  when  applied  to 
the  strip  and  the  cube.  Working 
equations  for  these  cases  are  repro- 
duced in  the  Appendix  to  this  paper. 
On  the  other  hand,  calculations  for 
the  long  cylinder  and  the  sphere  are 
little  more  difficult  than  the  one- 
dimensional  case.  The  finite  volume 
elements  may  be  taken  as  concentric 
cylindrical  or  spherical  shells,  respec- 
tively. Mass  balances  then  lead  to 
the  following  general  expressions : 

For  the  long  cylinder : 


AT, 


aSn"+1 


2a  +  lL 


(9) 

where  a  =  -y—  ,  rn  being  the  inside 

radius  of  element  n,  and  the  symbol 
Snm  signifies  the  integral: 

jTPm 

5V"=  /  Pip 

Jpn 

For  the  special  case  of  the  central 
core,  where  a  =  0,  Equation  (9) 
reduces  to: 

AO 


AT„  =  2SB 


s(&ry 


(10) 


and  for  the  outside  shell,  assumed  to 
have  only  half  the  thickness  of  the 
others, 

2 


A7\  = 


[aSS 


a  +  H 

+  k(Ar)(a  +  ]/2)(pa-p1)] 


Ad 


s(Ar)2 
(11) 


-25 

* 

O 

k-20 


PIG.  8.  EFFECT  OF  VARIABLE 
PERMEABILITY  ON  DRYING  RATE 

Slab,  half-thickness  4V4,  Ax  =  1,  (  =  1, 
Js  =  2.5  X  10-*,  vapor  pressure  as  in 
Fig.  3.  Curve  A,  P  =  4.4  X  10"*. 
Curve  B,  P  as  in  Fig.  9.  Curve  O,  P 
as  in  Fig.  10 
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For  example,  if  the  cylinder  is  di- 
vided into  five  concentric  elements, 
the  outermost  having  only  half  of 
the  standard  thickness,  Ar,  the  work- 
ing equations  are  as  follows  : 

ATX  =  [1.882SX2 

'+  2.118  k(A 

AT 2  =  (0.858  Si 


+  1.142  S21) 
AT&  =  (0.800  Sa* 

+  1.200  S32) 
ATt  =  (0.667  SJ 


Ad 


+  1.333  V) 
AT6  =  2(SB«) 


j(Ar)2 

AO 

s(Ar)2 

Ad 

s(Ar)2 

A6 


s(Ar)2 

The  mean  moisture  content  of  this 
cylinder  may  be  computed  as  follows : 

T  =  -^-(4.25  T,  +  7T2 

+  5r8  4-  3T4  4-  r„) 

For  the  general  element  of  the 
sphere : 

AT  -  a2Sn'+1+  (a+iyS^ 
a*  +  a+iA 

■  w  m 

For  the  central  element: 
ATm  =  35V-1  • 


Ad 


s(Ar)' 


(13) 

For  the  surface  element  (half- 
thickness)  : 

A7\  = 

„f ^S12+k(Ar)(a  +  y2y(pa-p1)-] 


a2  +  a/2  +  1/12 

Ad 
s(Ar)2 


(14) 


/ 


0  005  0.10  015 

MEAN  MOISTURE  CONTENT, T 


The  working  equations  for  a  sphere 
divided  into  five  concentric  elements 
are  as  follows : 


ATi  =  [1.770  Si2 
4-2.240ife(A 
AT2  =  (0.730  S2* 


+  1.300  Sz1)  J£r$ 


ATS  =  (0.632  Ss* 

+  1.420  Ss1) 
AT*  =  (0.429  SS 


j(Ar)s 

Ad 
s(Ar)' 


+ 1.714  st*) 

'  s(Ar)2 
Ad 


1 


s(Ar)'< 


T  =  -g^(27.2T1  +  37T2  +  19T3 
+  7Ti  +  T6). 

Precision  of  Approximate  Method 

Partial  tests  of  the  precision  of  the 
numerical  method  of  computation 
have  been  made  by  comparison  with 
the  results  of  analytical  solution  in 
certain  simple  cases.  Tables  1  and  2 
show  that  the  accuracy  of  estimation 
of  mean  moisture  content  is  amply 
good  for  most  purposes. 

Experience  with  the  routine  com- 
putation of  a  considerable  number 
of  cases  indicates  that  an  experienced 
calculating  machine  operator  can 
complete  about  fifty  steps  per  day  in 
a  problem  involving  only  one-dimen- 
sional diffusion,  the  solid  body  being 
divided  into  five-volume  elements, 
and  the  permeability  and  vapor-pres- 
sure functions  being  given  as  em- 
pirically plotted  curves.  Simple  com- 
putations such  as  those  leading  to 
Tables  1  and  2  proceed  considerably 
faster. 

Effect  of  Diffusion  Variables 
on  Drying  Rate 

It  may  not  be  apparent  immedi- 
ately that  the  drying  phenomena  that 
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FIG.   9.     FIRST  ARBITRARILY  ASSUMED 
PERMEABILITY 

would  be  described  by  Equation  (1) 
are  not  identical  with  those  corre- 
sponding with  Equation  (2)  if  the 
diffusivity  be  taken  as  constant  in  the 
one  case  and  the  permeability  as  con- 
stant in  the  other.  Figures  1  and  2 
show  successive  moisture  gradients 
through  the  half-section  of  an  infinite 
slab  for  the  two  cases.  Figure  3 
shows  the  vapor-pressure  isotherm 
that  was  used  in  calculating  the 
gradients  for  Figure  2 ;  it  is  that  of 
one  sample  of  dehydrated  carrots  at 
60°  C.  Figure  4  shows  the  relation 
of  drying  rate  to  mean  moisture  con- 
tent for  the  two  cases,  and  Figure  5 
the  relation  between  mean  moisture 
content  and  the  moisture  content  at 
the  center  of  the  slab.  The  sigmoid 
shape  of  some  of  the  gradients  in 
Figure  2  is  striking.  Gradients  having 
this  shape  have  been  found  experi- 
mentally by  several  workers  (5) 
(9)  (16)  (19)  (26)  (31)  (34) 
but  Ceaglske  and  Hougen  believed 
the  phenomenon  to  be  "inconsistent 
with  the  behavior  of  diffusion  in 
slabs  (9)."  Evidently  it  is  consistent, 
however,  with  a  diffusion  equation 
that  is  not  too  rigidly  restricted  by 
simplifying  assumptions. 

Variation  of  Diffusivity 
and  Permeability 

The  measurements  of  moisture 
transfer  through  cellophane  reported 
by  Doty,  Aiken,  and  Mark  (10)  do 
not  lend  themselves  to  precise  calcu- 
lation of  permeability  as  a  function 
of  moisture  content  because  the  hu- 
midity gradients  used  were  so  great 
that  a  wide  range  of  moisture  con- 
tents must  have  been  present  within 
a  single  sheet.  The  relation  shown 
in  Figure  6  may  be  deduced,  however, 
as  at  least  consistent  with  the  data. 
Figure  7  shows  the  corresponding  re- 
lation between  diffusivity  and  mois- 
ture content,  calculated  from  Figure 
6  and  a  vapor-pressure  isotherm  for 
viscose  rayon  reproduced  in  Interna- 
tional Critical  Tables.  (17)  The  state- 
ment made  earlier  that  the  choice  of 
potential  and  conductance  terms  is 


5 

O 
x 


□  I 

CD 

< 

£  0.5 

e 
ui 

O. 

0.2  1  A  

0         25         50         75  100 

VAPOR  PRESSURE,  mm.  Hg. 

FIG.  10.  SECOND  ARBITRARILY  ASSUMED 
PERMEABILITY 

mainly  a  matter  of  convenience  may 
be  illustrated  from  these  data.  All 
solutions  of  the  following  form  of 
the  diffusion  equation : 


will  be  identical  with  the  solutions  of 
Equation  (3),  namely: 

il  =  ±  *LCP*L\ 

dd       s       dx\     dx )' 

if  the  values  of  D  from  Figure  7  are 
used  in  the  former  and  the  values  of 
P  from  Figure  6  in  the  latter. 

The  effects  of  the  following  vari- 
ables on  drying  rate  in  the  lower  part 
of  the  falling-rate  phase  are  studied 
readily  by  the  use  of  the  approxi- 
mate numerical  procedure. 

Change  of  permeability  with  mois- 
ture content. 

Form  of  the  vapor-pressure  iso- 
therm. 

Change  of  the  ratio  between  sur- 
face film  coefficient  and  permeability. 

Change  in  the  shape  of  the  body 
(practically  limited  to  the  slab,  the 
cylinder,  and  the  sphere) . 

Change  in  the  thickness  of  the 
body. 

Figures  1,  2,  4,  and  5,  already  dis- 
cussed, show  the  variations  that  may 
be  expected  to  follow  from  the  use 
of  different  characteristics  for  the 
diffusion  potential  term  if  an  un- 
varying conductance  term  be  as- 
sumed. • 

Change  of  Permeability  with 
Moisture  Content 

Figure  8  compares  computed  dry- 
ing rates  for  the  infinite  slab,  assum- 
ing three  quite  different  permeability 
functions.  All  are  based  on  vapor 
pressures  as  the  potential  term,  and 
the  isotherm  is  that  reproduced  in  Fig- 
ure 3.  For  curve  A  the  permeability 
was  assumed  to  remain  unvarying 
with  moisture  content,  while  for 
curves  B  and  C  the  functions  repre- 
sented in  Figures  9  and  10,  respec- 
tively, were  assumed.  The  order  of 


magnitude  of  the  changes  was  sug- 
gested by  Doty,  Aiken,  and  Mark's 
results  for  cellophane,  but  the  func- 
tions were  given  simple  geometrical 
shapes  so  that  the  relation  between 
the  shape  of  these  curves  and  that  of 
the  derived  rate  curves  could  be  eas- 
ily analyzed. 

Figure  1 1  shows  a  group  of  mois- 
ture gradients  in  the  slab  for  the 
conditions  of  curve  B  in  Figure  8. 
Gradients  like  that  at  time  100  have 
been  observed  experimentally  by 
many  investigators.  The  material  is 
"case-hardened." 

The  peculiar  shapes  of  the  rate 
curves  indicate  how  profoundly  the 
drying  rate  may  be  affected  by 
changes  in  permeability  that  are,  in 
general,  known  to  occur.  A  note  of 
caution  on  the  interpretation  of  ex- 
perimental drying  rate  curves  is 
strongly  suggested.  The  existence  of 
a  substantially  straight  segment  of 
the  falling-rate  curve  is  no  guarantee 
that  the  straight  line  may  be  extrapo- 
lated all  the  way  to  the  base  line. 
Again,  since  an  estimation  of  drying 
time  almost  invariably  is  the  ultimate 
object,  any  line  of  best  fit  which  is 
to  be  drawn  through  a  set  of  points 
representing  successive  experimental 
drying  rates  should  deviate  from  the 
points  by  a  smaller  and  smaller  aver- 
age distance  as  zero  rate  is 
approached. 

Rate  curves  B  and  C  in  Figure  8 
have  no  simple  relationship  to  the  per- 
meability functions  on  which  they 
were  based.  This  is,  of  course,  only 
what  might  have  been  expected;  at 
any  given  moment  the  drying  rate  is 
the  result  of  a  complex  internal  situa- 
tion which  finds  some  parts  of  the 
solid  still  in  the  high-permeability 
range  while  others  have  already  dried 
to  a  condition  of  low  permeability. 
The  effect  of  the  sharp  break  in  per- 
meability shown  in  Figure  9  at  a 
vapor  pressure  of  about  70  mm.  (cor- 
responding to  a  moisture  content 
T  =  0.0800)  is  apparent  in  the  rate 
curve  at  a  mean  moisture  content  far 
above  T  =  0.0800.  No  presently 
known  procedure  can  be  used  to 
carry  the  calculation  in  the  opposite 
direction — that  is,  to  deduce  quanti- 
tatively a  permeability  function  from 
a  given  drying  rate  curve. 

Form  of  the  Vapor-Pressure 
Isotherm 

Figure  12  compares  the  computed 
drying  rates  for  an  infinite  slab,  using 
three  different  vapor-pressure  iso- 
therms, the  permeability  function  be- 
ing assumed  as  in  Figure  9  in  all 
three  cases.  The  three  isotherms, 
which  are  drawn  in  Figure  13,  were 
taken  from  experimental  data  on 
three  different  samples  of  dehydrated 
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vegetables ;  the  temperature  is  as- 
sumed in  all  cases  to  be  140  de- 
grees F. 

Effect  of  Surface 
Film  Coefficient 

Figures  14  and  15  show  drying 
rates  for  an  infinite  slab  as  affected 
by  changes  in  the  ratio  between  the 
surface  film  coefficient  and  the  per- 
meability. In  Figure  14  the  perme- 
ability is  assumed  to  remain  un- 
changed at  all  moisture  contents, 
while  in  Figure  15  it  is  assumed  to 
change  twenty-fold  within  a  narrow 
range,  as  in  Figure  9..  By  compari- 
son of  the  two  figures  it  may  be  seen 
that  when  the  permeability  is  variable 
it  is  the  ratio  of  the  surface  coeffi- 
cient to  the  low  range  of  permeability 
values  that  is  controlling,  even  in  the 
higher  levels  of  moisture  content. 

Drying  Rates  of  the  Slab, 

Cylinder,  and  Sphere 

Figure  16  compares  drying  rates 
for  an  infinite  slab,  an  infinitely  long 
cylinder,  and  a  sphere,  all  of  the  same 
half-thickness.  The  vapor-pressure 
isotherm  of  Figure  3  and  the  per- 
meability function  of  Figure  9  were 
used  in  all  three  cases. 

The  classical  diffusion  theory,  as- 
suming a  constant  diffusivity,  pre- 
dicts that  in  the  lower,  essentially 
linear,  portions  of  the  rate  curves  the 
ratio  of  the  drying  rates  for  three 
such  bodies  will  be  approximately 
1:2.3:4.4.  The  introduction  of  a 
variable  permeability  destroys  the 
linear  character  of  the  curves  and 
leads  to  somewhat  erratic  values  of 
these  ratios,  depending  to  some  ex- 
tent upon  the  form  of  the  permeabil- 
ity and  vapor-pressure  functions.  In 
this  example  the  ratios  in  the  lower 
portions  of  the  curves  are  approxi- 
mately 1:2.7:4.9. 

Effect  of  Thickness 

Classical  drying  theory  predicts 
that  if  the  surface  film  resistance  is 
negligible  the  drying  rate  will  vary 
inversely  as  the  square  of  the  thick- 
ness. Figure  17  compares  the  dry- 
ing rates  of  two  slabs  whose  thick- 
nesses are  in  the  ratio  of  4i^:9j/2; 
the  surface  film  resistance  is  negli- 
gible (coefficient  infinite)  and  the 
permeability  remains  constant  at  all 
moisture  levels  (the  diffusivity  there- 
fore changes  markedly  with  moisture 
content).  Comparison  of  the  ordi- 
nates  of  the  two  curves  shows  a  con- 
stant ratio  of  4:46,  which  is,  within 
the  limit  of  error,  the  square  of  the 
ratio  of  thicknesses.  The  classical 
rule  is  therefore  valid  for  the  more 
general  case.  This  conclusion  may 
be  demonstrated  also  analytically. 

If  the  drying  rate  that  is  being 


measured  is  that  of  an  assemblage 
of  pieces  which  vary  somewhat  in 
thickness  (as  is  always  the  case  in 
the  drying  of  prepared  vegetables) 
the  resulting  rate  curve  is,  of  course, 
a  composite.  The  effect  will  be  to 
mask  such  sharp  changes  in  curvature 
as  appear  in  curve  B  of  Figure  8. 
For  example,  suppose  50  per  cent  by 
weight  of  a  mixture  consists  of  pieces 
which  alone  would  have  the  drying 
rate  of  curve  B,  30  per  cent  consists 
of  pieces  having  only  half  the  stand- 
ard thickness,  and  20  per  cent  con- 
sists of  pieces  having  only  one-quar- 
ter the  standard  thickness.  Figure  18 
shows  how  the  composite  drying 
rate  compares  with  the  standard  dry- 
ing rate,  Curve  B.  (This  calculation 
ignores  the  mutual  effects  of  different 
pieces,  and  therefore  can  only  hold 
for  a  very  lightly  loaded  tray.) 

Shapes  of  Typical 
Experimental  Curves 

Figure  19  shows  several  experi- 
mental rate  curves,  taken  at  random 
from  the  large  number  determined 
during  the  past  few  years  in  this 
laboratory.  In  these  experiments  the 
materials  were  dried  in  shallow  layers 
on  screen-bottom  trays,  using  cross- 
flow  of  air.  The  results  have  more 
general  significance,  however,  because 
in  this  low-moisture  region  the  man- 


ner of  exposure  to  the  air  stream  has 
little  or  no  effect  on  the  drying  rate. 
The  literature  of  drying  theory  pro- 
vides few  data  suitable  for  compari- 
son. Some  of  the  drying  experiments 
(on  riced  white  potatoes),  reported 
by  Marshall  (24)  were  carried  to  a 
moisture  content  below  10  per  cent, 
and  the  few  points  in  this  low-mois- 
ture region  indicate  curvature  of  the 
rate  curves  generally  similar  to  the 
curvature  appearing  in  Figure  19. 

The  shapes  of  these  curves  are 
strongly  reminiscent  of  some  of  the 
theoretical  curves  —  for  example, 
curve  C  in  Figure  8.  Variation  of 
permeability  no  greater  than  has  been 
demonstrated  already  for  cellophane 
and  horn  membranes  would  account 
without  difficulty  for  the  experimen- 
tal rate  curves. 

Conclusions 

Approximate  numerical  solutions 
of  the  general  diffusion  equation  for 
the  infinite  slab,  the  long  cylinder,  and 
the  sphere  have  a  useful  degree  of 
precision  without  being  excessively 
laborious. 

Reasonable  assumptions  about  the 
variation  of  permeability  lead  to  cal- 
culated moisture  gradients  and  dry- 
ing-rate curves  similar  to  those  found 
experimentally  by  many  investi- 
gators. 


FIG.  12.  EFFECT  OF  SHAPE 
OF  VAPOR-PRESSURE  ISO- 
THERM  ON   DRYING  RATE 

Slab,  half -thickness  4%, 
Az  =  l,  s  =  l,  k  infinite, 
permeability  as  in  Fig.  9. 
Curves  A,  B,  and  O  cor- 
respond to  isotherms  A, 
B,  and  C  in  Fig.  13 
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FIG.   18.     VARIATIONS   IN   SHAPE  OF 
VAPOR-PRESSURE  ISOTHERM 

Curves  A  and  B,  dehydrated  carrots, 
curve  0,   dehydrated   potatoes.  Tem- 
perature, 140  degrees  F. 

The  shapes  of  the  drying-rate 
curves  bear  no  simple  relation  to  the 
quantitative  variation  of  permeabil- 
ity, and  no  procedure  for  deducing 
the  latter  from  the  former  is  known. 

When  the  permeability  is  variable, 
the  effect  of  surface  film  transfer  on 
drying  rate  is  controlled  by  the  ratio 
of  the  surface  coefficient  to  the  low 
range  of  permeabilities. 

When  the  permeability  is  variable, 
the  ratios  of  the  drying  rates  of  a 
slab,  a  long  cylinder,  and  a  sphere, 
of  equal  thickness,  are  no  longer 
fixed,  but  depend  upon  the  form  of 
the  permeability  function. 

The  drying  rate  is  inversely  pro- 
portional to  the  square  of  the  thick- 
ness of  the  piece  if  the  surface  film 
resistance  is  negligible,  regardless  of 
the  form  of  the  permeability  func- 
tion. 
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FIG.     14.  -   EFFECT    OF  SURFACE 
FILM    COEFFICIENT,  PERMEABIL- 
ITY CONSTANT 

Slab,  half-thickness  4%,  Ax  =  l,  *  =1, 
P  —  4  A  X  10-*,  vapor  pressure  as  In 
Fig.  8.  Surface  film  coefficients,  top  to 
bottom,  infinity,  5  X  10-*,  2.5  X  10"«, 
1  X  10-« 


FIG.    15.     EFFECT   OF  SURFACE 
FILM  COEFFICIENT,  PERMEABIL- 
ITY VARIABLE 

Slab,  half -thickness  4tt,  Ax  =  l,  (  =  1, 
P  as  in  Fig.  0,  vapor  pressure  as  In 
Fig.  8.  Solid  curve,  k  infinite;  dashed 
curve,  k  =  2.5  X  10"* 


0  OAS  0.10  0.IS 


MEAN  MOISTURE  CONTENT ,  T 


FIG.  16.     EFFECT  OF  SHAPE  OF 
BODY  ON  DRYING  BATE 

Vapor  pressure  as  in  Fig.  8,  permeabil- 
ity as  in  Fig.  9,  *  =  1,  Ax  =  &r  =  1, 
k    infinite.    Slab    half-thickness  and 
cylinder  and  sphere  radii,  4% 
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FIG.  18.    EFFECT  OF  MIXTURE  OF  SIZES 
ON  DRYING  RATE 

Slab,  Az  =  1,  i  =  l,  P  as  in  Fig.  9, 
k  =  2.5  X  10-4,  vapor  pressures  as  in 
Fig.  8.  Curve  B,  half-thickness  4%. 
Curve  A,  mixture  (by  weight)  60  per 
cent  standard  thickness,  30  per  cent  V4 
standard,  20  per  cent  V*  standard 
Rate  assumed  to  vary  as  3/2  power  of 
thickness 


MEAN  MOISTURE  CONTENT  ,T 


were  made  by  A.  H.  Brown,  P.  W. 
Kilpatrick,  and  their  coworkers. 
Vapor-pressure  isotherms  were  deter- 
mined by  B.  Makower  and  Fred  Tal- 
burt.  The  author  is  indebted  also  to 
W.  D.  Ramage,  Fred  Stitt,  H.  K. 
Burr,  and  M.  E.  Lazar  for  numerous 
helpful  suggestions. 


Nomenclature 

(Any  consistent  units) 


a  = 


ratio  of  inside  shell  ra- 
dius to  radial  incre- 


ment = 


Ar 


D  = 

*  = 

P  = 
Pa  = 

Pn  = 


concentration  of  mois- 
ture, mass  per  unit 
volume 

diffusivity  (concentra- 
tion as  potential) 

surface  mass  transfer  co- 
efficient 

moisture  vapor  pressure 

partial  pressure  of  mois- 
ture vapor  in  ambient 
air 

vapor  pressure  at  ele- 
ment no.  n 

change  in  vapor  pressure 
in  an  increment  of 
time 


P  = 


Ar  = 
s  = 

C  m  — 


T  = 


TB  = 


Tn  = 
AT  = 

A*  = 

A0  = 
dT/dd  = 


permeability  (vapor 
pressure  as  potential) 

mean  value  of  P  over  in- 
terval m  to  n 

inside  radius  of  cylindri- 
cal or  spherical  shell 
no.  n 

increment  of  radius 

mass  of  dry  solids  per 
unit  volume 

abbreviation  for 

JfPm 
'  PdP 

moisture  content,  pounds 
of  water  per  pound 
dry  solids 

moisture  content  at  cen- 
ter line  of  body 

moisture  content  in  equi- 
librium with  ambient 
air 

moisture  content  at  ele- 
ment no.  n 

change  in  moisture  con- 
tent in  an  increment 
of  time 

increment  of  distance, 
measured  normal  to 
the  drying  face 

increment  of  time 

drying  rate 
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FIG.  19.     TYPICAL  EXPERIMEN- 
TAL    DRYING     RATES,  DEHY- 
DRATED VEGETABLES 

Curve  A,  Run  747,  carrot  half-dice, 
160  degrees  F.,  wet-bulb  90  degrees; 
Curve  B,  Run  397,  6/32-in.  potato 
strips,  160  degrees  F.,  wet-bulb  90  de- 
grees; Curve  O,  Runs  667-9,  %-in. 
potato  dice,  150  degrees  F.,  wet -bulb, 
90  degrees 


In  all  of  the  calculations  described 
in  this  paper  a  uniform  initial  distri- 
bution of  the  moisture  was  assumed, 
(T  =  0.20).  Any  other  desired  dis- 
tribution, for  example  parabolic, 
could  have  been  used  without  chang- 
ing the  procedure  or  increasing  the 
labor.  The  effect  of  the  initial  distri- 
bution vanishes  rapidly,  however,  and 
is  negligible  in  the  drying-rate  curves 
presented  in  this  paper. 
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Appendix 

Working  Equations  for  Approximate 
Calculation  of  Diffusion  in  the 
Square  Strip  and  the  Cube 

The  equations  are  derived  by  solv- 
ing for  AT  the  appropriate  mass  bal- 
ances for  each  finite  element  of  the 
body. 


PIG.    20.      TRIANGLE  REPRESENTATIVE 
OF  INFINITE  SQUARE  STRIP 


Infinitely  Long  Square  Strip 

It  is  evident  from  considerations 
of  symmetry  that  the  solution  is  need- 
ed only  for  the  right  triangle  shown 
in  Figure  20,  constituting  one-eighth 
cf  the  cross-section  area.  Figure  21 
shows  how  the  triangle  will  be  di- 
vided into  elements  if  the  surface 
element  has  only  half  the  thickness 
of  the  interior  ones,  and  indicates 
a  convenient  numbering  system  for 
the  elements.  No  diffusion  will  occur 
either  across  the  diagonal  side  or 
across  the  vertical  side. 

In  general,  the  mass  balance  for 
an  element  will  contain  four  terms, 
two  representing  flow  into  the  ele- 
ment and  two  flow  out  of  it.  The 
equation  for  a  typical  interior  ele- 
ment, say  No.  24,  will  be : 

AT24  =  (S2434  +  S2425  +  S2423 

Ad 

+  S24U)  s(Ax)2 

The  equation  for  a  typical  surface 
element,  say  No.  14,  must  take  ac- 
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count  both  of  the  smaller  volume  of 
the  element  and  the  smaller  diffusion 
area  from  two  of  the  adjacent  ele- 
ments : 

Ar14  =  2[s^  +  y2s^  +  y2s^* 

Ad  . 


+  k(Ax)(pa-p14)] 


s(Ax)' 


The  remaining  equations  may  be 
written  in  the  same  manner.  The 
mean  moisture  content  of  the  strip 
may  be  calculated  by  weighting  the 
moisture  content  of  each  element 
proportionally  to  the  volume  of  that 
element. 


The  Cube  . 

Figure  22  shows  the  right  tetra- 
hedron which  sufficiently  represents 
the  diffusional  process  in  the  cube, 
and  Figure  23  indicates  a  convenient 
numbering  system  for  the  elements. 
The  tetrahedron  has  only  1/48  of 
the  volume  of  the  entire  cube.  The 
pyramidal  element  at  the  corner, 
Number  111,  has  only  1/48  of  the 
volume  of  a  full  cubical  element  like 
Number  245,  while  the  larger  pyra- 
midal elements  along  the  diagonal, 
such  as  Number  333,  have  1/6  of 
the  volume  of  the  cubical  element. 

In  general,  the  mass  balance  for  an 
element  contains  six  terms,  three  rep- 
resenting inflow  and  three  represent- 
ing outflow.  For  an  interior  element, 
say  Number  234,  the  equation  will 
take  the  following  form : 

Ar234  =  (S2342"  +  5-234224  +  SW33 
.884 


+  .SW34) 


Ad 


s(Ax)' 


The  equation  for  a  surface  ele- 
ment, say  Number  133,  taking  ac- 
count of  the  change  in  the  element 
volume  and  of  differences  in  diffu- 
sional area,  is  as  follows : 

AT133  =  4[H5,i83238  +^5-1S31S4 

+  i/iW23  +  H*(M(A, 
A6 


-  P133)] 


s(Axy 


DIVISION  OF  TRIANGLE  INTO 
ELEMENTS 


The  laborious  nature  of  the  calcula- 
tions for  the  strip  and  the  cube  is 
obvious.  For  a  given  fineness  of  sub- 
division there  are  three  times  as  many 
elements  to  be  calculated  in  the  square 
strip  as  in  the  slab,  and  seven  times 
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FIG.    22.      TETRAHEDRON  REPRESENTA- 
TIVE OP  CUBE 


as  many  in  the  cube  as  in  the  slab. 
In  addition,  the  complexity  of  the  cal- 
culation is  greater  at  each  element. 
For  the  fineness  of  subdivision  shown 
in  Figures  21  and  23  a  good  machine 
computer  can  complete  only  about 
twelve  to  fifteen  steps  per  working 
day  for  the  square  strip,  and  only 
five  or  six  for  the  cube.  The  shapes 
of  the  drying-rate  curves  are  quite 
similar  to  those  of  the  long  cylinder 
and  the  sphere,  respectively. 

Presented  at  the  Regional  Meeting,  San 
Francisco,  California,  August  25-28,  1946. 


Discussion 

G.  N.  Harcourt  (Buffalo,  N.  Y.) : 
I  had  the  pleasure  of  reading  Mr. 
Van  Arsdel's  well-conceived  paper, 
and  I  was  very  much  interested  in  it. 

Unfortunately,  the  little  I  ever 
knew  about  differential  equations  has 
long  since  been  forgotten  so  that 
part  of  the  paper  has  been  accepted 
on  faith. 

This  kind  of  paper  is  of  great 
potential  value  since  it  guides  the 
thinking  of  others  who  have  specific 


problems  to  solve.  The  specific  prob- 
lem may  not  be  susceptible  to  solu- 
tion by  the  mathematical  processes 
described  in  this  paper  but  the  ideas 
that  Mr.  Van  Arsdel,  his  co-workers, 
and  other  workers  have  set  forth 
are  of  great  value. 

Our  engineers  have  known  for  a 
long  time  that  the  rate  of  removal 
of  moisture  from  material  containing 
high  percentages  of  moisture,  70  per 
cent  to  10  per  cent,  is  much  greater 
than  it  is  from  materials  containing 
only  10  per  cent  to  4  per  cent  mois- 
ture. For  example,  the  condenser  on 
a  vacuum  dryer  must  be  designed  to 
take  care  of  the  large  quantity  of 
vapor  which  is  given  off  by  the  wet 
material  at  the  beginning  of  the  cycle. 


PIG.    23.      DIVISION    OF  TETRAHEDRON 
INTO  ELEMENTS 


Table  3. 

Test                          A  B  C  D  E  F  G 

RPM  drums    6-14  6-'/2         10        13-^         11  9  10 

Steam  pressure  ...  44  60          60  60  60  60  60 

Lb./sq.ft./hr   4.01  4.59       5.77         9.33        6.62  5.31  5.74 

Per  cent  moisture 

in  dry  material  ..  2.62  0.94       1.395       10.39         2.90  1.43  1.78 
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It  is  therefore  our  practice,  when 
making  recommendations  to  our  cus- 
tomers, to  keep  these  possibilities  in 
mind  and  we  will  recommend  on  oc- 
casion the  installation  of  two  differ- 
ent kinds  of  dryers,  one  which  will 
remove  moisture  rapidly  down  to  the 
7-10  per  cent  point,  the  other  which 
will  provide  a  longer  retention  time, 
the  application  of  heat  in  a  different 
way  or  which  is  in  some  other  way 
more  suitable  for  the  final  reduction 
of  the  moisture  content,  rather  than 
to  attempt  to  do  the  job  on  one  type 
of  dryer. 

Much  of  our  experience  has  been 
with  drum  dryers  in  which  capacities 
are  expressed  in  pounds  of  dry  ma- 
terial per  square  foot  of  drum 
surface. 

At  one  time  we  dried  a  solution  of 
a  sodium  salt  under  the  following 
conditions  (see  Table  3  also). 

Moisture  in  solution  at  start  =  70 
per  cent. 

Steam  pressure  in  drums  =  60  per 
cent. 

Our  customer  had  specified  that 
the  dry  material  should  contain  a 
maximum  of  0.1  per  cent  moisture 
which,  as  the  foregoing  indicates,  we 
could  not  obtain  economically  by 
means  of  the  drum  dryer.  We  sug- 
gested finally  that  the  drum  dryer 
be  used  to  reduce  the  moisture  to  2-3 
per  cent  at  which  moisture  content 
the  material  from  the  drum  dryer 
could  be  handled  easily  in  a  conveyor 
and  did  not  cake.  This  partially  dried 
material  was  then  fed  to  a  direct- 
fired  rotary  hot-air  dryer  where  it 
was  subjected  to  a  very  high  tem- 
perature and  the  last  traces  of  mois- 
ture driven  off.  The  combination  of 
the  drum  dryer  and  rotary  dryer 
gave  an  economical  and  compact  in- 
stallation which  was  eminently  satis- 
factory to  the  user. 
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5.02 

5.82 

6.77 

6.26 

1.82  2.0        4.13-5.27  2.64-3.30 
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Has  any  consideration  ever  been 
given  to  the  effect  of  concentration  of 
the  soluble  components  in  the  juices 
in  vegetables,  fruits  and  other  cel- 
lular materials  while  being  dried.  At 
the  beginning  of  drying,  when  the 
material  contains  70-80  per  cent 
moisture,  the  solids  content  of  the 
juice  may  be  4-6  per  cent  but  as  dry- 
ing progresses  the  solids  content  will 
increase.  It  is  a  well-known  fact  that, 
generally,  as  the  solids  content  of  a 
solution  increases,  the  vapor  pres- 
sure of  the  water  vapor  from  the 
solution  will  decrease.  If  the  drying 
temperature  remains  constant,  the 
rate  of  removal  of  moisture  will 
decrease. 

R.  L.  Huntington  (University  of 
Oklahoma,  Norman)  :  I  notice  a  sim- 
ilarity between  the  moisture  and 
temperature  gradients  for  batch  cool- 
ing of  solids.  I  wonder  if  you  made 
a  correlation  between  heat  flow  and 
drying  rates. 

W.  B.  Van  Arsdel:  No,  we  have 
not.  The  special  point  we  had  in 
mind  was  the  effect  of  a  very  large 
change  in  permeability  and  so  far  as 
I  know  that  condition  is  not  likely 
to  be  met  in  the  study  of  heat  flow. 
I  doubt  that  temperature  gradients 
in  a  body  undergoing  cooling  would 
ever  exhibit  such  peculiar  shapes  as 
the  moisture  gradients  in  a  body 
undergoing  drying.  Conductivities 
and  heat  capacities  are  more  nearly 
independent  of  temperature. 

R.  L.  Huntington:  You  didn't 
make  any  attempt  to  study  the  effec- 
tive permeability  of  the  flow  of  vapor 
compared  to  that  of  liquid  ? 

W.  B.  Van  Arsdel:  We  begged 
that  question.  It  will  require  care- 
fully planned  experiments  to  deter- 
mine whether  the  diffusion  of  water 
under  these  conditions  is  actually 
proportional  to  a  concentration  dif- 
ference or  to  a  vapor-pressure  differ- 
ence. In  an  isothermal  system  either 
concept  can  be  used,  and  the  results 
will  be  identical.  We  chose  to  de- 
velop our  results  on  the  vapor-pres- 
sure basis,  because  we  believe  that 
will  be  more  useful  when  we  come  to 
consider  nonisothermal  systems. 

Written  Discussion 

G.  M.  Dusinberre  (Virginia 
Polytechnic  Institute,  Blacksburg)  : 
The  author  points  out  that  it  is  im- 


possible to  infer  from  experimental 
drying-rate  curves  the  details  of  the 
phenomena  which  combine  to  pro- 
duce the  observed  result.  It  is  pos- 
sible only  to  test  the  validity  of  var- 
ious hypotheses  by  comparing  their 
logical  conclusions  with  the  experi- 
mental facts.  The  author  has  not 
been  content  to  use  the  oversimplifiedi 
classical  analysis  but  has  taken  a 
worth-while  step  in  testing  more  com- 
plex hypotheses,  using  the  numerical 
analysis  as  a  tool. 

The  writer  has  a  philosophic  ob- 
jection to  the  use  of  the  word  "ap- 
proximate" in  the  title  of  this  paper 
and  elsewhere,  as  characterizing  the 
numerical  treatment  by  contrast  with 
the  differential  treatment.  Engineers 
and  scientists  use  differential  equa- 
tions so  much,  and  to  such  great  ad- 
vantage, that  there  is  a  tendency  to 
join  the  mathematicians  in  regarding 
these  equations  as  fundamental,  rath- 
er than  the  physical  situation  which 
the  equations  are  set  up  to  represent. 

The  present  physical  problem  is  to 
predict  the  effects,  observable  on  the 
large  scale,  of  discrete  molecules 
moving  in  a  heterogeneous  and  dis- 
continuous region.  As  one  approxi- 
mation we  pass  over  the  discon- 
tinuities, treat  the  region  as  homo- 
geneous, and  apply  the  rules  of  cal- 
culus. As  another  approximation  we 
set  up  arbitrarily  discontinuities  at 
convenience  and  apply  the  rules  of 
arithmetic.  Both  methods  are  ap- 
proximate, physically.  Both  are  ex- 
act, mathematically,  in  the  sense  that 
there  is  a  demonstrable  relation  be- 
tween the  initial  assumptions  and  the 
final  solution.  But  the  assumptions 
are  slightly  different  in  the  two  cases. 
It  is  only  when  the  one  solution  is 
compared  with  the  other  set  of  as- 
sumptions that  it  deserves  to  be  called 
an  approximation. 

No  one,  least  of  all  the  writer,  can 
overlook  the  great  advantage  of  an 
analytical  solution,  when  one  is  ob- 
tainable. This  field  of  mathematics 
needs  no  apology  of  defense.  Neither 
does  it  need  the  undeserved  halo  of 
being  called  "exact"  in  its  physical 
application.  The  numerical  treat- 
ment evidently,  at  this  time,  needs 
some  defense  and  explanation.  But 
it  does  not  deserve  to  be  character- 
ized as  merely  approximate,  by  con- 
trast with  the  other. 

For  the  engineer  and  scientist  the 
test  of  accuracy  is  agreement  with 
experiment.  From  this  viewpoint, 
the  author  has  done  some  excellent 


work  in  demonstrating  the  accuracy., 
of  the  method  which  he  refers  to  still 
as  approximate. 

Incidentally,  it  may  be  borne  in 
mind  that  an  analytical  solution  is 
usually  in  the  form  of  an  infinite 
series.  If  the  first  few  steps  of  a 
numerical  solution  are  inaccurate,  so 
are  the  first  few  terms  of  a  series.  If 
the  numerical  solution  may  run  into 
a  divergence,  so  may  a  series.  But 
there  are  always  remedies. 

The  author  stated  that  the  writer's 
methods  were  unsatisfactory  when 
variation  of  diffusivity  was  wide  and 
abrupt.  This  offers  difficulty  with 
any  treatment,  but  it  may  be  pointed 
out  that  .the  JPdp  used  by  the  author 
is  exactly  analogous  to  the  jkdT  used 
by  the  writer.  Thus,  whatever  ad- 
vantages or  disadvantages  are  found 
must  lie  in  the  details  of  application, 
since  the  two  treatments  are  the  same 
in  principle. 

W.  B.  Van  Arsdel:  Professor 
Dusinberre's  point  about  my  designa- 
tion of  the  numerical  method  as  "ap- 
proximate" is  very  well  taken.  The 
approximations  with  which  we 
should  be  seriously  concerned  are 
those  concealed  within  the  simplify- 
ing assumptions.  I  suspect  that  many 
of  us  are  unduly  impressed  by  the 
"elegance"  of  analytical  solutions. 
Numerical  solutions  look,  in  contrast, 
clumsy  and  unimaginative;  uncon- 
sciously, then,  we  assume  they  must 
be  less  precise  representations  of  the 
true  physical  situation. 

The  comment  made  in  the  paper 
about  Professor  Dusinberre's  semi- 
graphical  method  for  handling  va- 
riable-diff  usivity  calculations  was 
perhaps  not  sufficiently  explained. 
The  principle  is,  as  he  says,  exactly 
the  same  as  the  one  I  have  applied. 
However,  if  a  radical  change  in  dif- 
fusivity or  permeability  is  involved 
(such  as  the  20-  to  40-fold  changes 
discussed  in  the  paper)  the  construc- 
tion of  his  "correction  chart"  to  a 
sufficient  degree  of  precision  be- 
comes exceedingly  laborious  in  itself. 
More  pertinent  from  the  standpoint 
of  the  exploratory  work  reported 
here  is  the  fact  that  many  different 
combinations  of  permeability  func- 
tion and  vapor-pressure  isotherm 
were  being  investigated.  Each  dif- 
ferent combination  would  have  re- 
quired a  new  correction  chart.  If  a 
large  number  of  calculations  were  to 
be  made  for  a  single  substance  the 
labor  invested  in  preparation  of  the 
chart  might  be  justified. 
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